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. - , Abstract 

o- ■ 

, Two approaches to the study of cosmological density perturbations in modified theories of Palatini gravity 

' have recently been discussed. These utilise, respectively, a generalisation of Birkhoff 's theorem and a 

\ direct linearization of the gravitational field equations. In this paper these approaches are compared 

I . and contrasted. The general form of the gravitational lagrangian for which the two frameworks yield 

^ ' identical results in the long- wavelength limit is derived. This class of models includes the case where the 

, lagrangian is a power-law of the Ricci curvature scalar. The evolution of density perturbations in theories 

l/^ ■ of the type f{R) = R — c/R^ is investigated numerically. It is found that the results obtained by the 

' two methods are in good agreement on sufficiently large scales when the values of the parameters {b, c) 

, are consistent with current observational constraints. However, this agreement becomes progressively 

^ ' poorer for models that differ significantly from the standard concordance model and as smaller scales are 

considered. 

u ; 

1 Introduction 

CN ■ There is now overwhelming observational evidence that the universe is currently undergoing a phase of 

. accelerated expansion. This evidence arises from high redshift supernovae surveys [TJ [2J [31 El [5] , observations 

' of large-scale structure O [7], baryon acoustic oscillations ^, and high-precision data from the Cosmic 

^ ! Microwave Background (CMB) pfTUj. 

' An accelerating universe poses a great challenge to modern cosmology. R is very difficult to explain 

ly-j . such behaviour within a conventional general relativistic framework. The simplest way to generate a phase 

' of acceleration is to introduce a cosmological constant into the field equations. Although this is entirely 

. consistent with the available data, the microphysical origin of such a term remains a mystery. 

' Ri view of this, a large number of alternative models have been proposed (see for a recent review). 

^ . In most cases, these can be classified into two broad groups: those that invoke an exotic matter source 

for the dark energy [12] and those that modify the gravitational sector of the theory [13]. Examples of 
the latter include generalised theories of gravity based on non-linear functions, f{R), of the Ricci scalar R. 
' Such modifications to the (linear) Einstein-Hilbert action typically arise in effective actions derived from 

string/M-theory Tl [JH [H [IT] . 

Given the large number of candidates within both approaches, it is important to determine whether 
observations can be employed to discriminate between them. In general, the homogeneous dynamics encoded 
in the modified Friedmann equations of f{R) theories can always be expressed in terms of a conventional 
relativistic cosmology sourced by an effective perfect fiuid. Thus, it is necessary to include inhomogeneities in 
order to lift the degeneracy between the two frameworks. This can be achieved, for example, by considering 
the evolution of density perturbations and a study of perturbation theory in /(i?) gravity is therefore of 
considerable interest. Two approaches have been developed recently in this context. One possibility is 
to follow the standard procedure employed in relativistic perturbation theory and derive the perturbation 
equations by linearizing the gravitational field equations |18'ill9j. On the other hand, an alternative procedure 
has recently been put forward by Lue, Scoccimarro and Starkman (LuSS) |20| which employs a generalised 
version of Birkhoff's theorem (see also Ref. [21 ). This procedure has the benefit of greatly simplifying 
the analysis, but suffers from the drawback that the degree of its applicability in more general settings is 
presently not known in detail. 

The purpose of the present work is to perform a detailed comparative study of the evolution of pertur- 
bations obtained by employing the LuSS procedure and the direct linearization of the field equations. Such 
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a comparison can serve as a crucial step in clarifying the status of the LuSS approach in non-linear gravity 
theories. 

When studying generalised f{R) gravity, a choice must be made as to which independent fields should 
be varied in the action. There are two possibilities, which are referred to as the 'metric' and 'Palatini' 
frameworks, respectively. In the former case, only variations with respect to the metric are considered, 
whereas the action is varied with respect to both the metric and the connection in the Palatini framework 
|22j . Both approaches result in identical field equations for the Einstein-Hilbcrt action, where f{R) (x R. 
More generally, however, the metric approach results in fourth-order equations, whereas the Palatini variation 
generates a second-order system [23]. Moreover, the theories of the type f{R) = R — c/R^ based on metric 
formalism have difficulty in passing the solar system tests [24] (see, however, [25]) and producing the correct 
Newtonian limit [55]. Such theories also suffer from gravitational instabilities as discussed in [57]. Along 
with these concerns, a recent study has found that these types of theories can not produce a standard matter 
epoch |28| . In this paper, we consider theories based on the Palatini variational method which can avoid the 
problems outlined above. Such theories have received considerable attention in recent years as candidates 
for explaining the present-day acceleration of the universe [HI [22] |23l |29l |30l [Ml |32l |33l |34] . 

The outline of the paper is as follows. In section 2, we briefly review the Palatini variational method and 
present the perturbation evolution equations derived by using the LuSS procedure and the direct linearization 
of the field equations. In section 3, we derive the necessary conditions on the general form of the gravitational 
lagrangian f{R) for the two approaches to be compatible. We perform a numerical analysis in section 4 to 
quantitatively compare the two frameworks for the class of f{R) theories where the action includes an inverse 
power of the Ricci curvature scalar. We conclude with a discussion in Section 5. 



2 f(R) gravity in the Palatini formalism 
2.1 The Field Equations 

We consider the class of non-linear gravity theories defined by the action 



^ - 7^ / d^x^ [/(i?)] + Sra{g, ^), (1) 
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where k = SttG is a constant, f{R) is a differentiable function of the Ricci scalar R = g^^'^ R^^{g,V) and 
Rf^^{g, F) is the Ricci tensor of the affine connection P^^. The matter action Sm is a functional only of the 
metric tensor and matter fields, ip. 

In the Palatini framework, the affine connection and the metric are treated as independent variables. 
Extremizing action ([T]) with respect to the metric tensor yields the gravitational field equations (see, e.g., 
Ref. [12]): 

F{R)R^, - ^/(i?)5M- = ^^T^., (2) 

where FiR) = df/dR and 



V = -^77£ (3) 



2 SS„ 
^Sgt^" 

defines the energy-momentum tensor. Contraction of Eq. ^ yields the algebraic constraint equation 

RF{R) - 2f{R) - kT, (4) 

where T = T^^. 

On the other hand, varying action ([T]) with respect to the affine connection P^^ and contracting implies 
that [55] 

V,[FiR)^gn = 0, (5) 

where V is the covariant derivative defined by P^-y. The solution to Eq. ^ is given by writing P^^ as the 
Levi-Civita connection for a new metric /i^t^ = F{R)g^i,^, which is conformally equivalent to the spacetime 
metric g^n,. As a result, we may write the Ricci tensor of the affine connection in the form 

i?p.(f ) = R,A9) + ^ - j^^uF - _5^,V.VV, (6) 
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where R^j^^ig) is the Ricci tensor of the Levi-Civita connection. It can then be shown that the field equations 
([2]) and constraint equation ^ can be derived from a scalar-tensor action of the form [35] 
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(7) 



where the scalar field (f) = F and V{(t)) = R{4))F — represents the potential. Action ([7]) corresponds 

to the Brans-Dicke theory with a dilaton-graviton coupling loq = —3/2. 

We assume throughout that the background spacetime is given by the spatially flat and isotropic Friedmann- 
Robertson- Walker (FRW) metric 

ds^ = a{T)^{-dT^ + S^jdx'dx^), (8) 

where a(r) is the scale factor of the universe and r = J dt/a{t) defines conformal time. We further assume 
that the matter content of the universe is represented by a pressureless perfect fluid with energy-momentum 
tensor T^^, — diag(— p„i, 0, 0, 0). The Ricci tensor ([6]) may then be employed to derive the generalised 
Friedmann equation for this cosmology: 



6F H 




(9) 



where a dot denotes differentiation with respect to coordinate time t and H = d/a deflnes the Hubble 
expansion parameter. After substituting the trace equation ([4]), the generalised Friedmann equation can be 
rewritten in the form 



3f-RF 
6F 



1 - 



3 F^ni2f^RF) Y 
2 F{F - RF^r) J 



(10) 



where a comma denotes differentiation. 



2.2 The Perturbation Evolution Equations 

In conventional cosmology, there exists an interesting equivalence between the Newtonian and general rel- 
ativistic frameworks. Both approaches result in identical background evolution equations (i.e. Friedmann 
equations) as well as evolution equations for the scalar perturbations. The former coincidence results from 
the fact that there is an analogue of Newton's sphere theorem in general relativistic settings, i.e., Birkhoff's 
theorem holds. The correspondence for the perturbation evolution arises in the absence of vector and tensor 
fluctuations. 

Recently, a procedure has been put forward by Lue, Scoccimarro and Starkman [20] which relies on the 
assumption that this Newtonian analogy, including Birkhoff's theorem, holds in the more general setting of 
modified gravity theories. According to this procedure, it is assumed that the growth of large-scale structure 
can be modelled in terms of a uniform sphere of dust of constant mass, such that the evolution inside the 
sphere is determined by the FRW metric. Using Birkhoff's theorem, the spacetime metric in the empty 
exterior is then taken to be Schwarzschild-like. The components of the exterior metric are then uniquely 
determined by smoothly matching the interior and exterior regions. 

The overdensity 6(t) of the spherical distribution of pressureless matter with mass M and radius r is 
defined by 

l + 6{t)^—^, (11) 

where p{t) represents the background energy density. The matching conditions imply that f = r{H^ + H) 
and the evolution of the density perturbation is then given by i36j 




S + 2HS - \2H + — ] S = 0, (12) 

or, equivalcntly, by 

6" + HS' - (^-2H']S^0, (13) 
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where a prime denotes difTerentiation with respect to conformal time and 7Y = aH = a. Eq. (I13p can also 
be derived by assuming that the continuity and Friedmann equations apply directly to the fluctuations 21J. 

Recently, the evolution of perturbations in f{R) gravity was investigated using the LuSS procedure 
pm [3T] . The advantage of this approach is that the growth of the density contrast can be expressed in terms 
of a single quadrature involving the Hubble parameter and the scale factor j20[ [37] : 



S.HjJL^^ (14) 

In principle, therefore, the evolution of the perturbations can be determined once the background dynamics 
has been specified. 

However, the validity of the LuSS procedure has yet to be established in generalised gravity. It is 
important, therefore, to compare this approach with the method that directly linearises the gravitational field 
equations. The corresponding evolution equation for comoving, linear density perturbations in a pressureless 
universe was recently derived by Koivisto and Kurki-Suonio (KKS) using this direct method and found to 
have the form |TB] 

^ + 3FnH2Fn + F') ^ 

GF^n^jn" - 2n'n) + GF'^njn^ ~ W) + F'F{;m"n ^ m'^ - n'^k^) + QF"Fn{n' ~ n^) 

3FH^{2Fn + F') ~ ' 

where k is the comoving wavenumber that arises due to the Fourier decomposition. 

We will refer to Eqs. and (fT5|) as the LuSS and KKS perturbation equations, respectively. We will be 
interested in identifying the domain where the equation based on the LuSS procedure provides an accurate 
description for the evolution of the perturbations. In the following section, we adopt an analytical approach 
with the aim of identifying the general form of the gravitational lagrangian, f{R), for this to be the case. 



3 Analytical Results 

A direct comparison between the LuSS equation (fT3|) and the KKS equation (fT5|) suggests that the latter 
should be rewritten in the form 

d" + ^ns'-c(j^-2n'y^o, (16) 

where the parameters ^ and ( are defined by 

2FF"n - 2F'^n - 2FF'n' 

^^^^ Fn^{2Fn + F') ' ^ ' 



^^'^^ U" - 2H'H^^ ^^^^ 3{2FH + F'){H" - 2H'n) ^ ' ""^^^ 



and 



respectively. The form of Eq. (|T6l) implies that the LuSS and KKS equations are equivalent when ~ C, = \, 
but it is clear that this occurs only for Einstein gravity where F' = 0. Indeed, the most striking difference 
is the presence of the gradient term in the KKS equation. Such a term also arises in the corresponding 
density perturbation equation derived in the metric variational approach [38 . The origin of this term can be 
understood from the dynamical equivalence between Palatini gravity and Brans-Dicke theory, as expressed 
in Eq. ([7|). Fluctuations in the pressureless matter induce perturbations in the scalar field (i.e. the 
Ricci curvature), which in turn generate a pressure gradient in the fluid. In general, the sound speed of the 
fluctuations in the cold dark matter is given by 

2 

" 3{2Fn + F') ■ ^^^^ 

The magnitude of ^ is independent of k and is therefore unaffected by the specific choice of scale. 
However, C, contains a gradient term which is proportional to fc^ and this may be significant on small scales. 
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Consequently, the evolution of the perturbations will indeed be different in the two approaches. However, 
the gradient term becomes negligible in the long- wavelength limit (which corresponds formally to — > 0). 
In this limit, a necessary and sufficient condition for equivalence between the LuSS and KKS equations is 
that ^ = 1 and this constraint is satisfied when 

FF"n - F'^n - FF'n' = o. (20) 

Eq. (|20p may be viewed as a second-order, non-linear differential equation for F(t). One solution to this 
equation is that of general relativity with a cosmological constant, f{R) = R — A. More generally, if F' 7^ 
and F" 7^ 0, we may define a parameter Y = F' /F. This reduces Eq. (pO| to the remarkably simple form 

Y' H' 

— = — , (21) 
which admits the integral Y ~ YqH., where Yq is an arbitrary integration constant. This in turn implies that 

F^Foa^", (22) 

where Fq is a second integration constant. 

On the other hand, the trace equation ([4]) for a universe sourced by pressureless matter reduces to the 
condition [31] 



ao.[2f-R-^) . (23) 

Hence, substitution of Eq. into Eq. yields a first-order, non-linear differential equation in the 

gravitational lagrangian f(R): 

T^VYs/ - ^l;') = constant, (24) 



dRj V dR^ 
where n = S/Yq. 

Eq. (I24p is a particular example of d'Alembert's equation and may be solved in full generality [55] . 
Since we are interested in the functional dependence of the lagrangian on the Ricci scalar, we may rescale 
/ without loss of generality such that the constant on the right-hand side of Eq. is unity. If we now 
define the functions 

M^lf, iV.ir#V" (25) 



2dR 2 \dR^ 

and denote p = df /dR, Eq. ([M)) can be expressed in the form f{R) — RM{p) + N{p). Differentiating this 
expression with respect to R then yields 



. ^ dp 
P = M(p) + - 



^ dM{p) ^ dNjp) 
dp dp 



(26) 



However, Eq. (|26p can be expressed as a linear differential equation in the dependent variable R and 
independent variable p: 

^ _^ ^ !!_ (27) 

dp p p2+n- ) 

Hence, solving Eq. ([?7|) by the method of integrating factors yields the general solution to Eq. ([M)) in a 
parametric form: 

f-l^P+^n^ (29) 
where Cq is an arbitrary integration constant and P is a free parameter. 
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Eqs. P51) - (P5)1 represent the general form of the gravitational lagrangian f{R) for the LuSS and KKS 
equations to be compatible in the long wavelength limit. It is interesting that for this class of theories the 
sound speed of the fluctuations is constant with a numerical value given by 

When Co = 0, which is equivalent to the asymptotic limit where R is sufficiently small, the gravitational 
action depends on a simple power of the Ricci scalar: 

/(i?) cx (31) 

For this class of theories the Friedmann equation PH]) reduces to 



o 3 + 2n 3 



-2 



which in turn implies that the background dynamics is given by a power-law solution for the scale factor, 
a oc 7Y-2"/(3+n) ^ ^2n/(3+n) |;^qj ^ Consequently, the cosmic dynamics is equivalent to that of a conventional 
relativistic universe dominated by a perfect fluid with a constant equation of state. Finally, the parameter 
C simplifies in this case to 

^ 3(l + n)(3 + n)(3 + 2n)H2' ^ ' 

In conclusion, therefore, the above analysis indicates that the LuSS equation should provide a good 
approximation to the full evolution equation for the linear density perturbation on sufficiently large scales in 
any modified gravity theory that asymptotes in the low-energy limit to a power-law in the Ricci curvature 
scalar. On the other hand, for fixed values of n and 7i, the LuSS equation becomes progressively less accurate 
as we move to smaller scales (i.e. as k increases). In the following section, we will quantify these conclusions 
further by performing numerical calculations for a specific class of modified gravity theories. 



4 Numerical Results 

Motivated by the results of the previous Section, we consider the class of gravity theories defined by 

f{R) = R-^, (34) 

where b and c are free parameters whose values are constrained by observations. Such theories have recently 
been considered as possible candidates for explaining the late-time acceleration of the universe [32l [33l [34] . In 
particular, a recent study found that data obtained from CMB, baryon oscillation and large-scale structure 
observations constrains the parameters (6, c) to lie in the ranges h E [—0.2,1.2] and c e [—3.5,6.6] at the 
68% confidence level [31]. The best-fit model corresponds to the values (6, c) = (0.027, 4.63) and the ACDM 
concordance model is represented by (&, c) = (0, 4.38). These values are consistent with the results of other 
studies that employ CMB and supernovae data [31]. 

For the above choice of parameters, we have made a detailed comparative study of the evolution of the 
density perturbations for both the LuSS equation and the KKS equation The results of such a 

comparison can be quantified by defining a 'fractional difference' parameter 

. SluSS — 5kKS .or:\ 

A= 7 , (35) 

okks 

where subscripts 'LuSS' and 'KKS' refer to the results obtained using the LuSS and KKS equations, re- 
spectively. Thus, the two approaches are completely compatible when A = 0. This parameter is defined in 
such a way that the difference between the two approaches is of the same order as the KKS approach when 
A ~ 0{1). To a first approximation, therefore, it is reasonable to suppose that the LuSS equation becomes 
unreliable when A sa 1. 
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There are three physical parameters in the field equations whose values need to be specified in the 
numerical integrations. These are ^mo, Rq and Hq^ where a subscript zero indicates present-day values 
and ri„i is the normalised matter energy densitj0. However, only two constraint equations are available, 
corresponding to the Friedmann equation (jlOp and the trace equation (jj]). In order to be consistent, therefore, 
we specify the value of Hq to be unity, as is the usual practice (see, e.g., [31]). We then use the constraint 
equations ^ and pU|) to determine flmo and Rq. The choice of Eq. (|35p implies that the initial value of 
the perturbation 6 is unimportant. Finally, we need to specify the scale of the perturbations. By fixing 
the wavenumber at a particular value, one focuses on perturbations that entered the horizon at a particular 
epoch. For illustrative purposes we consider the values k = 5 and k = 20, corresponding to scales which 
remain within the horizon throughout our numerical evolution. 

The left hand panel of Fig. [1] illustrates the evolution of A when c — 4.38 and fc = 5, with b taking values 
in the range b € [0, 1]. As expected, A = for the ACDM concordance model (given by 6 = 0), since it is 
known that the LuSS equation is exact in this case. On the other hand, increasing the value of b causes the 
behaviour of the two approaches to deviate and the quantitative difference becomes more pronounced as b 
is increased. 




Figure 1: Illustrating how the fractional difference parameter. A, varies with scale factor a as 6 is increased. 
Here c — 4.38 and values of b are assigned to each curve. The case & = corresponds to the ACDM model. 
The left hand panel corresponds to fc = 5 whereas the right hand panel corresponds to A: = 20. 

We have verified that our results remain qualitatively similar when the parameter values lie in the ranges 
b € [—0.2, 1.2] and c G [—3.5, 6.6], respectively. An important outcome of these results is that for values of 
the parameters consistent with recent observations, the agreement between the LuSS and KKS approaches 
is good in the sense that A < 0.1 for 6 < 0.2. This implies that the LuSS equation provides a good 
approximation to the full (linear) perturbation theory (for this value of k). This can be understood by 
noting that observations constrain theoretical models to lie close to the ACDM point, where it is known that 
the LuSS equation is exact. 

Further inspection of the left hand panel of Fig. [1] indicates that as the value of b is increased, the models 
take longer to move away from the ACDM point A = 0, but those with smaller values of b subsequently find 
it easier to approach A = at later times. We may gain further insight into the origin of this behaviour by 
investigating the evolution of the quantity Q = 1 — F. This vanishes at all times for Einstein gravity but is 
given by Q = —bcR~^^^''^ for the class of models (p4| . This parameter therefore provides a measure of the 
deviation away from general relativity. Our numerical calculations indicate that initially R w 0(10'^) and, 
consequently for larger values of b, the scale factor must grow to a larger value before the Ricci scalar has 
fallen sufficiently for the correction term Q to become dynamically significant. In other words, the onset 
of acceleration occurs at later times for larger b. On the other hand, the correction term in f{R) that is 
proportional to R~^ will become more important as the universe expands. The analysis of Section 3 then 
indicates that the accuracy of the LuSS equation will improve as f{R) asymptotes to a power-law form. 

^Notc that in modified gravity theories of the type considered here, this parameter need not necessarily be unity in a spatially 
flat universe. 
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Consequently, A will begin to decrease back to zero at later times. 

We find qualitatively similar behaviour at larger values of k. The right hand panel of Fig. [1] illustrates 
the corresponding evolution of A when k = 20. As expected, models with lower values of b move away from 
the A = point at smaller values of the scale factor. The model with the lowest non-zero value of 5 = 0.2 
crosses the solutions for b — 0.4 and b = 0.6. This can be understood from Eq. ([33| . which implies that the 
magnitude of C, depends on the ratio k^/H^ — k'^/cP. At a formal level, therefore, increasing the value of k 
is equivalent to ending the numerical calculation at a fixed k but with a smaller value for the scale factor. 

However, the quantitative agreement between the solutions of the LuSS and KKS equations is poor when 
k = 20 and A rapidly exceeds unity in this case. This discrepancy arises primarily because the deviation of 
the parameter C. away from unity is more pronounced at larger k. Fig. [2] illustrates the evolution of ( for 
the different values of k. 
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Figure 2: Illustrating the evolution of the parameter C defined by Eq. p8)) in the text for the parameter 
values k — h (left panel) and /c = 20 (right panel). The LuSS procedure for the evolution of the perturbations 
becomes progressively less accurate as the deviation of this quantity from unity becomes more pronounced. 



5 Conclusion 

In this paper, we have studied the evolution of density perturbations in generalised theories of gravity where 
the field equations are derived via the Palatini variational approach. We focused on models where the 
energy-momentum tensor is sourced by a pressureless perfect fluid. Two approaches to the study of density 
perturbations have recently been developed in the literature [T51 I?!]. These involve, respectively, an 
application of Birkoff's theorem to modified gravity (the LuSS method) and the linearization of the full 
field equations (the KKS approach). In the former case, the evolution of the perturbations is determined 
entirely by the background dynamics and no pressure gradients are present in the perturbation evolution 
equation. However, such terms do arise in the linearization approach, which takes into account the fact that 
perturbations in the fluid induce fluctuations in the Ricci curvature which in turn modify the sound speed 
of the fluctuations in the matter. 

In the long-wavelength limit, these gradient terms are negligible. We have identified the most general /(i?) 
theory of gravity, as summarised in Eqs. (pS)) and (P^)) . for the LuSS and KKS approaches to be compatible 
in this limit. A particular case of this class of theories arises when /(i?) is a simple power law of the Ricci 
curvature scalar. This is interesting because such terms are expected to arise generically as corrections to 
the Einstein-Hilbert action at low energies. Furthermore, theories of this type result in a background scaling 
solution, in the sense that the homogeneous dynamics is equivalent to that of a conventional relativistic 
cosmology where the pressure and energy density of the perfect fluid redshift at the same rate. It would 
be interesting to explore whether this scaling behaviour is a necessary condition for compatibility between 
the LuSS and linearization methods in more general theories of modifled gravity. For example, a power-law 
cosmology arises in the Palatini variation of Ricci squared gravity, where / oc (i?^''i?^j^)"/^ [IT] . 
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We numerically investigated a specific class of power-law theories of the type (|34p and compared the 
LuSS and KKS approaches on smaller scales where gradient terms become significant. We found that when 
the parameters of the underlying theory take values that are consistent with cosmological observations, the 
LuSS procedure provides a reasonably good approximation to the complete linearised theory if k is not too 
large (i.e. of the order of a few or less). However, the agreement between the two approaches soon breaks 
down on smaller scales. 
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